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Figure 


ILLUSTRATIONS 


The  total  tracking  force  with  the  electric  and  magnetic 
components,  calculated  for  an  electron  beam  using  the 
analytic  model,  are  plotted  versus  the  beam  coordinate, 
t.  The  10-kA  beam  with  a  0.25-cm  displacement  propa¬ 
gates  through  a  1-cm  radius  channel  with  0.2-cm_l  con¬ 
ductivity. 

The  tracking  force,  on  a  10-kA  beam  with  a  15-cm  rise 
time  in  a  1-cm  channel  with  a  =  0.5  cm~l,  is  plotted 
using  the  results  from  the  analytic  model  and  IPROP. 

In  (a),  the  total  force  is  given  and,  in  (b),  the  elec¬ 
tric  and  magnetic  components  of  the  force  are  plotted. 
Note  the  good  agreement  when  the  forces  are  relatively 
large. 

The  tracking  force  on  a  10-kA  beam  with  a  0.25-cm  dis¬ 
placement  is  calculated  using  the  analytic  model  with 
X  =  1,  2  and  4  and  compared  with  Lee’s  expression.  The 
1-cm  channel  has  a  -  1  cm“^.  The  drift-tube  radii  are 
2.7,  7.4  and  55,  corresponding  to  X  =  1,  2  and  4,  respec¬ 
tively. 

The  maximum  tracking  force  on  a  10-kA  beam  in  a  1-cm 
channel  with  a  =  0.2  cm"!  is  plotted  for  variable  beam 
rise  time. 
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The  tracking  force,  calculated  using  IPROP,  is  sketched 
versus  t.  The  channels  considered  are  conductivity, 
density  and  combination  (conductivity/density).  The 
10-kA,  50-MeV  beam  has  a  0.5-cm  radius  with  a  Bennett 
profile.  The  Gaussian-shaped  channel  has  r^  =  1  cm. 

The  electric  and  magnetic  components  of  the  total 
tracking  force  are  plotted  for  a  10-kA  beam  with  a 
15-ns  rise  time  in  a  1-cm  conductivity  channel. 

The  electric  and  magnetic  components  of  the  total 
tracking  force  are  plotted  for  a  10-kA  beam  with  a 
15-ns  rise  time  in  a  1-cm  density  channel. 

The  electric  and  magnetic  components  of  the  total 
tracking  force  are  plotted  for  a  10-kA  beam  with  a 
120-ns  rise  time  in  a  1-cm  density  channel. 

The  tracking  behavior  of  a  trumpet  beam  and  a  2-cm 
radius  pencil  beam  is  compared.  The  3-cm  conductivity 
channel  was  offset  0.75  cm.  The  trumpet  shape  is 
given  by  Eq.  33. 
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ILLUSTRATIONS  (Concluded) 


The  tracking  forces  for  a  trumpet  beam  in  a  density 
channel  is  plotted.  The  3-cm  radius  channel  had  a 
0.1-atm  on-axis  density. 
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Optimization  of  channel  tracking  in  conductivity  channel 
Optimization  of  tracking  in  density  channel 
Optimization  of  tracking  in  density/conductivity  channel 
Trumpet  beam  tracking  in  conductivity  channel 


CHAPTER  1 
INTRODUCTION 


In  this  paper,  we  will  study  the  response  of  an  intense  relativistic 
electron  beam  to  a  preformed  channel.  The  channel,  which  may  have  initial 
conductivity,  depressed  density  or  both,  could  be  produced  by  an  initial 
charged-particle  beam,  laser  or  electrical  discharge.  If  a  beam  is  intro¬ 
duced  displaced  from  the  channel  axis  and  through  its  own  interaction  with 
the  channel  is  forced  towards  the  axis,  the  beam  is  said  to  "track". 
Maxwell's  complete  equations  will  be  used  in  analytic  and  numerical  work  to 
study  channel  tracking  in  detail. 

The  new  analytic  work  presented  in  this  paper  builds  upon  the  electro¬ 
static  tracking  model  of  Lee.1  Lee  first  developed  an  analytic  tracking 
model  for  a  beam  propagating  entirely  inside  a  square  conductivity  chan¬ 
nel.  He  showed  that,  in  the  limit  of  small  conductivity,  the  problem 
becomes  purely  electrostatic.  In  order  to  take  into  account  the  magnetic 
effects  associated  with  more  sizable  conductivity,  we  used  the  full  frozen- 
field  equations  to  derive  the  tracking  force  for  this  simplified  problem. 
The  analytic  expression  illustrates  the  effects  of  both  electrostatic  and 
magnetic  tracking.  As  the  conductivity  approaches  zero,  Lee's  expression 
is  obtained. 

The  numerical  simulations  were  carried  out  by  the  three-dimensional 
simulation  code  IPROP.2  IPROP  is  a  particle-in-cell  code  with  air  chemis¬ 
try  based  on  the  Lawrence  Livermore  National  Laboratory  algorithm3  as  modi¬ 
fied  by  Science  Applications  International  Corporation  (SAIC).4  The  code 
was  used  in  the  ultrarelativistic  limit  for  a  beam  which  is  not  allowed  to 
respond  dynamically  (slug  beam).  In  this  case,  the  relatively  small  forces 
associated  wit'1  tracking  were  easily  identified.  IPROP  employs  the  full 
Maxwell  equations  without  approximations  such  as  those  suggested  by  Lee5 
and  agrees  well  with  the  derived  analytic  result.  The  complete  equations 
are  necessary  for  the  detailed  study  of  channel  tracking,  since  there  is  a 
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delicate  balance  between  electrostatic  and  magnetic  forces  near  the  head  of 
the  beam. 


As  is  shown  later,  the  forces  involved  in  channel  tracking  are  only 
dominant  fairly  close  to  the  head  of  the  beam.  Both  magnetic  and  electro¬ 
static  components  of  the  tracking  force,  observed  in  the  numerical  simula¬ 
tions,  decay  away  several  centimeters  into  the  beam.  The  electrostatic 
force  modeled  by  Lee1  is  created  by  the  beam  potential  acting  on  the  chan¬ 
nel  electrons.  This  produces  a  dipole  force  which  tends  to  push  the  beam 
back  toward  the  channel  axis  when  the  beam  is  displaced.  The  magnetic 
tracking  force  is  due  to  the  radial  and  azimuthal  currents  of  the  channel 
electrons  which  are  driven  by  the  beam  potential  as  well  as  to  the  axial 
displacement  current  in  Maxwell's  equations.  Since  the  distance  from  the 
beam  head  that  a  beam  tracks  has  only  a  slight  dependence  on  the  rise  time, 
short  rise  times  and,  hence,  large  displacement  currents  are  desirable. 

For  these  beams  tracking  is  dominated  by  the  magnetic  forces.  Electro¬ 
static  forces  dominate  for  longer  rise  times.  It  is  for  shorter  rise  times 
( <4  ns)  that  the  full  field  equations  are  most  needed. 

Previous  numerical  work  by  Masamitsu6  and  by  Hui  and  Lampe7  using 
Lee's  ultrarelati vistic  equations5  has  given  somewhat  smaller  values  than 
IPROP  for  the  tracking  force  in  the  short  rise-time  cases.  IPROP  has  cal¬ 
culated  forces  an  order  of  magnitude  greater  for  a  0.5  ns  rise-time  beam 
than  calculated  by  Hui  and  2-3  times  greater  than  Masamitsu.  For  a  rise 
time  of  4  ns,  IPROP  agrees  reasonably  well  with  Hui's  results.8  Another 
interesting  difference  seen  in  the  numerical  simulations  for  a  conductivity 
channel  is  the  optimal  initial  conductivity  for  tracking.  IPROP  studies 
have  indicated  values  as  high  as  3  (where  conductivity  is  normalized  as 
4ira^a'/c,  where  ab  is  the  beam  Bennett  radius  and  o'  is  the  scaler  conduc¬ 
tivity).  Previous  results  have  suggested  a  conductivity  of  0.1.  The 
higher  values  for  conductivity  and  tracking  may  be  accounted  for  by  the 
different  modeling  of  the  magnetic  fields  in  IPROP. 
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In  Chapter  2,  we  will  present  analytic  work  which  includes  the  track 
ing  force  derivation  for  the  complete  frozen-field  equations.  The  numeri 
cal  computations  involving  IPROP  simulations  are  presented  in  Chapter  3. 
In  Chapter  4,  we  discuss  the  physical  interpretations  of  the  results  and 
give  our  conclusions  in  Chapter  5. 
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CHAPTER  2 
ANALYTICAL  MODELS 


The  electrostatic  and  magnetic  dipole  forces,  acting  on  a  displaced 
electron  beam  in  a  conductivity  channel,  nearly  cancel  close  to  the  head  of 
the  beam.  This  is  the  crucial  region  which  will  determine  whether  or  not 
the  entire  beam  will  track  the  channel,  since  the  strong  pinch  force 
generally  will  cause  the  rest  of  the  beam  to  follow  the  head.  The  complete 
field  equations  are  necessary  to  fully  understand  what  is  happening  in  this 
region.  In  this  section,  we  will  consider  a  sharp-edged  beam  which  is 
propagating  entirely  inside  a  square  conductivity  channel.  A  brief  deriva¬ 
tion  of  the  resulting  tracking  force  on  the  beam  is  given.  The  result  is 
then  used  to  check  the  accuracy  of  the  code  IPROP  and  compared  with  Lee's 
expression.  All  quantities  will  be  given  in  cgs  units  with  the  tracking 
force  given  occasionally  in  gauss. 

A.  FROZEN-FIELD  EQUATIONS 

Maxwell's  equations  in  an  axially  moving  frame,  when  separated  into 
forward-  and  backward-propagating  components,  are  given  by,9 
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where  c  and  v  are  the  speed  of  light  and  the  reference  frame  velocity, 
respectively,  and  a  =  4ira'/c. 

The  frozen-field  equations  are  obtained  by  setting  v  =  c  and  dropping 
time  derivatives.  With  some  manipulation,  the  frozen-field  equations 
become 

Er  =  f(lfEz  *f  kh)  ■  <7 
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with  t  =  -z/c. 

When  a  is  constant,  Eqs.  11  and  12  simplify  greatly  into  uncoupled 
second-order  equations: 


s7  77/  2 


Furthermore,  Eq.  13  can  be  integrated  yielding 


We  now  have  six  equations,  Eqs.  7-10,  14  and  15,  valid  for  a  constant 
conductivity  region  in  the  frozen-field  limit. 


B.  DERIVATION  OF  THE  TRACKING  FORCE 

We  will  solve  the  above  equations  for  a  channel  with  a  square  conduc¬ 
tivity  profile  (constant  a  inside  r  =  a  and  zero  elsewhere).  The  tracking 
force  is  obtained  for  an  electron  beam  propagating  totally  inside  the  chan¬ 
nel.  The  equations  are  Fourier  analyzed  (simply  replace  second-order  dif- 

2 

ferentials  with  respect  to  9  with  -m  )  in  order  to  determine  the  dipole 
force  acting  on  the  beam  (force  in  the  m  =  1  mode).  The  only  approximation 
will  be  to  neglect  the  ad/dt  term  in  Eqs.  14  and  15  for  r  <  a.  For 
o  a  ln(b/a)  <  4,  where  b  is  the  drift  tube  radius,  this  term  was  found  to 
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be  important  only  at  the  boundary  of  the  channel  and  must  be  kept  in  the 
boundary  conditions.  The  sudden  build-up  of  surface  charges  and  currents 
at  r  =  a  is  responsible  for  the  significance  of  the  term,  particularly  near 
the  beam  head.  To  reduce  the  error  in  calculating  the  dipole  force,  we 
will  analytically  subtract  out  the  forces  symmetric  about  the  beam  by 
separating  the  fields  into  beam-  and  channel -symmetric  components. 


We  can  quickly  solve  the  m  =  1  mode  equations  outside  the  channel  with 
no  approximation,  since  the  ad/dt  term  is  identically  0.  The  m  =  1  fields, 
where  the  superscript  (1)  denotes  m  =  1  quantities,  are 


Ez(1)  «  Bz(1)  =  a  y  (r/a)'1 
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for  b  »  a,  where  y  and  6  are  x-dependent  only. 

The  boundary  conditions  are  determined  from  the  complete  frozen-field 
equations  (Eqs.  7-12)  by  evaluating  them  at  r  =  a.  The  useful  boundary 
conditions  to  the  derivation  are 


Ez  continuous  at  a 


(17a) 


B  continuous  at  a 


(17b) 


which  results  from  the  continuity  of  E0^  -  BQ^.  The  advantages  of  the 


last  boundary  condition  will  be  more  obvious  later. 


Taking  od/dx  =  0  inside  the  channel,  we  can  immediately  determine  B_ 


from  Eq.  14, 


Bz(1)  =  a  ^  y  (r/a) 


where  the  continuity  of  Bz  has  been  used.  Splitting  the  fields  into  beam- 


and  channel -symmetric  parts,  we  express  Ez  as 


=  a  A^(r)  +  a  ~  a  (r/a)  , 


where  a  is  dependent  only  on  x.  This  is  helpful,  since  only  the  second 
term  in  Eq.  19  need  be  considered  when  calculating  the  tracking  force. 


To  obtain  A^  ,  we  first  consider  the  case  where  both  the  beam  and 
channel  are  symmetric  about  the  z  axis  (m  =  0).  From  Eq.  15,  the  solution 
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where  Iz  is  the  total  current.  The  m  =  1  solution  is  obtained  by  dis¬ 


placing  the  beam  by  e.  When  e  is  small,  we  get 
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We  can  now  solve  for  a  by  noting  6  drops  out  of  the  fourth  boundary 
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condition  (Eq.  17d)  and  then  solving  for  a  in  terms  of  the  known  Av  '  (a). 
After  a  great  deal  of  algebra,  we  obtain 
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X  =  o2a2  4n(b/a)  . 


(22b) 


The  tracking  force  for  a  small  beam  displacement  is  given  by 

Ft  =  7  I  lEr  -  Be  -  E9  -  Br)111  Jz  rdr  •  (23) 

where  only  terms  involving  a  need  to  be  considered.  Again,  after  much 
algebra,  we  obtain  the  tracking  force  per  current: 
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where  we  adopt  the  convention  that  a  positive  force  indicates  tracking. 

The  tracking  force  is  proportional  to  e/a  and  is  independent  of  the  other 
spatial  characteristics  of  the  beam.  Eq.  24  is  valid  up  to  about  X  =  4, 
where  the  od/dx  =  0  approximation  breaks  down. 

C.  RESULTS 

The  physics  of  channel  tracking  is  best  understood  by  separating  the 
tracking  force  into  electrostatic  and  magnetic  components.  These  are 
obtained  from  Eq.  22  and  are  given  by 
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The  electrostatic  force  behavior  is  described  by  the  build  up  of  charge  at 
the  channel  boundary  and  its  subsequent  decay.  In  Eq.  25,  the  first  term 
in  the  integral  controls  the  build  up  of  charge  at  the  channel  boundary. 
This  term  decays  as  exp(-or).  The  second  term  controls  the  charge  decay 
from  r  =  a,  and,  for  X  =  0,  decays  more  slowly  as  exp(-<rr/2).  Thus,  in  the 
limit  of  X  =  0,  the  charge  builds  up  indefinitely  for  positive  dl^/d-r. 

The  first  and  second  terms  in  the  magnetic  force  equation  correspond 
closely  to  those  of  Eq.  25.  The  first  term  describes  the  effects  of  the 
transverse  outward  currents,  which  are  responsible  for  the  growth  of  the 
channel  surface  charge,  and  the  dipole  axial  displacement  current 
(dE  ^  /dt)  caused  by  the  resulting  asymmetric  channel  electrons.  Both 
effects  are  responsible  for  the  transverse  magnetic  fields  which  cause 
tracking.  The  second  term  controls  the  transverse  return  currents  which 
drain  the  surface  charge  and  cause  detracking.  Electrostatic  and  magnetic 
forces  are,  therefore,  intimately  related.  Note  that,  for  X  >  1,  the  solu¬ 
tions  to  Eqs.  25  and  26  become  oscillatory.  In  this  case,  the  transverse 
currents  alternate  direction  as  they  decay. 

Now,  we  will  consider  a  1 0 -k A  beam  offset  1/4  cm  in  a  1-cm  channel 
with  constant  a  =  0.2  cm~*.  The  time-dependent  beam  current  is  given  by 
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where  1/6  is  the  current  rise  time  which  we  will  take  as  15  cm  here.  The 
resulting  tracking  forces,  F  ,  Fb  and  the  total  force,  are  plotted  in 
Figure  1  for  b  =  25  cm.  We  see  at  the  beam  head  (t  =  0),  the  magnetic 
force  is  dominant  with  the  electric  force  giving  detracking  (negative 
force).  Later  on,  as  the  channel  electrons  are  able  to  redistribute  them¬ 
selves  (this  takes  place  on  time  scale  1/a),  the  electrostatic  force  begins 
to  dominate.  As  expected,  this  occurs  at  t  =  5  cm  in  Figure  1.  The  mag¬ 
netic  force  decays  and  eventually  becomes  detracking  at  t  =  12  cm.  If  the 
conductivity  is  sufficient,  the  magnetic  detracking  will  eventually 
dominate. 

We  found  that  IPROP  (using  the  same  channel  geometry  and  beam  as  the 
analytic  model)  is  in  good  agreement  (within  10%)  with  the  derived  expres¬ 
sion  for  the  total  tracking  force  for  X  <  4.  In  Figures  2a  and  2b,  the 
tracking  force,  on  the  same  beam  as  above,  is  calculated  using  IPROP  and 
the  analytic  model.  In  this  case,  the  1-cm  radius  channel  has  a  0.5-cm”1 
conductivity  (X  =  1).  In  Figure  2a,  the  total  tracking  forces  calculated 
using  IPROP  and  the  model  are  very  close  particularly  near  the  beam  head. 
When  the  force  is  separated  into  its  electric  and  magnetic  components 
(Figure  2b),  the  agreement  is  good  until  the  magnetic  component  becomes 
relatively  small.  This  deviation  results  from  the  neglect  of  the  ad/dr 
terms  in  Eqs.  14  and  15.  This  term  may  be  important  to  an  individual  com¬ 
ponent  of  the  force  if  its  magnitude  is  small  compared  to  the  total  force. 
For  X  >  4,  neglecting  the  ad/dx  term  is  no  longer  valid,  and  the  total 
force  calculated  analytically  deviates  from  IPROP.  The  good  agreement  of 
IPROP  to  the  rigorously  derived  analytic  expression  for  smaller  X  gives  us 
confidence  in  the  accuracy  of  IPROP,  which  will  be  used  extensively  in  the 
next  section. 

We  will  now  compare  our  new  expression  with  Lee's  for  the  above 
example  with  varying  X  (actually,  the  conductivity  and  the  channel  radius 
will  remain  constant  while  we  vary  the  drift  tube  radius).  The  total 
tracking  force  obtained  from  Eq.  24  for  current  profile  27  is 
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Figure  1.  The  total  tracking  force  with  the  electric  and  magnetic 
components,  calculated  for  an  electron  beam  using  the 
analytic  model,  are  plotted  versus  the  beam  coordinate. 
The  10- kA  beam  with  a  0.25-cm  displacement  propa¬ 
gates  through  a  1-cm  radius  channel  with  0.2-cm-l 
conductivity. 
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Figure  2.  The  tracking  force,  on  a  10- kA  beam  with  a  15-cm  rise 
time  in  a  1-cm  channel  with  e  =  0.5  cm'1,  is  plotted 
using  the  results  from  the  analytic  model  and  IPROP. 
In  (a),  the  total  force  is  given  and,  in  (b) ,  the 
electric  and  magnetic  components  of  the  force  are 
plotted.  Note  the  good  agreement  when  the  forces 
are  relatively  large. 
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For  the  same  case,  Lee's  expression  is  given  by. 


r/"  -  (a /2)e" 
B  -  a/2 


It  should  be  noted  that  Eq.  28  limits  to  Lee's  expression  in  the  limit  of 
small  X.  This  is  illustrated  in  Figure  3  where  Eq.  28  is  plotted  for 
X  =  1,  2  and  5,  while  holding  a  =  1  cm'1  and  a  =  1  cm.  Since  Lee's  model 
neglects  magnetic  effects  which  become  important  for  increasing  X,  Eq.  28 
deviates  from  Lee's  expression  at  large  X.  At  early  times  (small  t),  the 
magnetic  tracking  of  the  ful 1 -equation  expression,  seen  in  Figure  1, 
results  in  larger  tracking.  Later,  as  the  dipole  charge  builds  and  the 
electrostatic  component  becomes  dominant,  the  total  force  approaches  that 
of  the  purely  electrostatic  expression  of  Lee.  Note  the  distinct  oscilla¬ 
tion  in  the  tracking  force  (for  X  =  5)  which  was  predicted  earlier  for 

X  >  1.  In  this  example,  we  varied  the  drift  tube  radius  to  change  X.  How¬ 

ever,  the  same  deviation  is  seen  when  o  is  increased  and  a  and  b  are  held 
constant. 

Finally,  we  consider  the  rise-time  dependence  of  Eq.  28.  In  Figure  A, 
we  see  the  peak  tracking  force  plotted  versus  rise  time.  As  the  rise  time 

decreases  to  30  cm,  the  tracking  force  rises  sharply.  However,  as  the  rise 

time  approaches  the  inverse  of  the  channel  conductivity  (5  cm),  the 
increase  in  tracking  is  slowed.  This  is  caused  by  the  channel  electrons 
being  unable  to  respond  to  the  sharp  current  rise.  This  creates  an  optimal 
force  of  76  gauss  as  the  rise  time  approaches  zero. 
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The  tracking  force  on  a  10- kA  beam  with  a  0.25-cm  dis¬ 
placement  is  calculated  using  the  analytic  model  with 
X  =  1,  2  and  5  and  compared  with  Lee's  expression.  Th 
1-cm  channel  has  -  =  1  cm_l.  The  drift-tube  radii  are 
2.7,  7.4  and  150,  corresponding  to  X  =  1 ,  2  and  5, 
respecti vely . 
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Figure  4 


The  maximum  tracking  force  on  a  1 0 - k A  beam  in  a  1-cm 
channel  with  -  =  0.2  cm-1  is  plotted  for  variable 
beam  rise  time. 


The  derived  analytic  expression  illustrates  the  balance  between  elec 
trostatic  and  magnetic  forces  as  well  as  giving  us  a  check  on  I  PROP ‘ s 
behavior.  The  key  addition  to  the  complete  frozen-field  equations  is  the 
axial  displacement  current  which  drives  large  dipole  magnetic  fields  near 
the  beam  head.  These  fields  guide  the  beam  back  towards  the  z  axis.  IPR 
agrees  well  with  the  analytic  model  using  the  same  beam/channel  geometry. 
We  anticipate,  therefore,  that  in  a  realistic  beam/channel  configuration 
where  magnetic  tracking  dominates,  IPROP  should  yield  reliable  results 
which  may  differ  markedly  from  those  of  Lee's  ultrarelativistic  equations 


CHAPTER  3 

NUMERICAL  SIMULATIONS 


The  results  from  the  code  I  PROP  studies,  involving  a  more  complex 
channel  and  electron  beam,  are  presented  in  this  section.  The  three  types 
of  channels  considered  are  conductivity,  density  and  conductivity/density 
(combination).  The  radial  profiles  of  the  channels  are  Gaussian  with  con¬ 
ductivity  and/or  density  given  in  terms  of  the  channel -symmetric  coordinate 
r'  as 

ct  =  aQ  exp(-(r'/rc)2)  (30a) 

and 

n  =  namb  -  nd  exp(-(r7rc)2)  ,  (30b) 


where  ag,  namb  and  n^  are  the  maximum  channel  conducti vity,  ambient  pres¬ 
sure  in  atmospheres  and  the  depth  of  the  density  channel,  respectively. 


The  beam  has  a  50-MeV  energy  (the  energy  is  used  only  in  the  air 
chemistry  routine  in  IPROP)  and  a  Bennett  profile  with  current  density 
given  by 


4  (  Z  Z  \ 
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where  Jg  is  the  on-axis  current  density.  We  will  consider  both  pencil-  and 
trumpet-shaped  beams.  The  r-dependent  current  profile  is  again  given  by 
Eq.  27.  The  beam  will  not  be  allowed  to  respond  dynamically  to  the  channel 
forces . 


The  simulation  code  IPROP  is  finite  differenced  in  the  z  and  r  di  rec 
tions  and  Fourier-analyzed  in  M  number  of  modes  in  the  e  direction.  We 


have  run  the  code  with  both  channel  and  beam  off  the  z  axis.  The  two 
cases  agree  to  within  10%  for  small  displacements  (<0.5  a^).  We  found  it 
advantageous  to  run  the  code  with  the  channel  offset,  since  a  smaller 
number  of  8  modes  and,  hence,  less  computer  time  and  memory  are  needed. 

The  number  of  modes  is  dependent  on  the  displacement  divided  by  the  width 
of  the  object  displaced  (the  channel  or  the  beam),  and  the  channel  is 
generally  wider  than  the  beam.  For  moderate  channel  offsets  (<1.5  r  ),  2-4 
modes  are  sufficient  to  determine  the  tracking  force  accurately  in  the 
channel -offset  case.  In  this  case,  the  tracking  force  was  also  less  sensi¬ 
tive  to  b  (drift  tube  radius  set  at  25  cm  here)  than  in  the  offset  beam 
case. 


The  conductivity  in  IPROP  evolves  due  to  impact  ionization,  electron 
avalanche  (generated  by  the  electric  fields)  and  recombination  based  on  the 
SAIC  algorithm.  An  energy-dependent  momentum-transfer  cross  section  also 
is  used.  A  critical  numerical  parameter  in  the  code  is  the  cell  size  in 
the  z  direction.  We  found  that,  for  the  initial  conductivities  used 
[<b  cm-1),  the  results  were  insensitive  to  change  in  the  cell  size  for 
values  less  than  0.25  cm. 

A.  PENCIL  BEAM 

We  first  consider  a  typical  pencil-beam  case  for  the  three  channel 
types.  The  beam,  again,  has  a  10-kA  peak  current  with  a  15-cm  rise  time 
and  ab  =  0.5  cm.  The  1-cm  channels  are  offset  3/4  cm  and  have  0.2  cm-1 
conductivity  and  n^  =  0.9  atm  (a  minimum  density  of  0.1  atm  for  n^b  =  1 
atm).  In  Figure  5,  the  tracking  force  on  a  beam  slice  (defined  in  Eq.  23) 
per  current  is  plotted  for  the  three  cases.  For  all  cases,  the  beam  tracks 
for  approximately  5  cm,  with  the  conductivity  channel  tracking  longest. 

The  highest  tracking  is  found  in  the  conductivity/density  channel  at  15 
gauss. 

The  density  channel  differs  from  the  conductivity  channel  in  that  the 
tracking  is  almost  totally  magnetic  with  the  significant  electric  forces 
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The  tracking  force,  calculated  using  IPROP,  is  sketched 
versus  :.  The  channels  considered  are  conductivity, 
density  and  combination  (conduc ti vi ty /dens i ty ) .  The 
1'0-kA,  50-MeV  beam  has  a  0.5-cm  radius  with  a  Bennett 
profile.  The  Gaussi an-shaped  channel  has  r  =  1  cm. 


being  detracking  through  6  cm.  Figures  6  and  7  break  the  tracking  force  of 
the  conductivity  and  density  channels,  respectively,  into  its  electrostatic 
and  magnetic  components.  Note  the  tracking  force  is  completely  magnetic  in 
the  conductivity  channel  up  to  4  cm  into  the  beam,  although  the  increasing 
electric  component  is  responsible  for  the  longer  tracking.  In  the  density 
channel,  the  electric  force  is  detracking  except  at  the  very  head  where  it 
reaches  only  0.4  gauss.  This  magnetic  dominance  is  typical  for  short  rise¬ 
time  beams. 

The  channel  electrons  have  time  to  respond  and  cause  electrostatic 
tracking  when  the  rise  time  increases.  Figure  8  sketches  the  tracking 
behavior  of  a  120-cm  rise-time  beam  for  a  conductivity  channel.  The  con¬ 
tribution  of  the  electric  forces  are  now  becoming  comparable  to  the  mag¬ 
netic  forces.  Both  the  electric  and  magnetic  components  cause  detracking 
at  12  cm  with  the  return-current  magnetic  detracking  dominating.  In  the 
density  channel  case,  the  magnetic  forces  will  generally  dominate  since  the 
conductivity  is  too  localized  near  regions  of  high  electric  field  to 
significantly  rearrange  the  channel  electrons  and  cause  electrostatic 
tracking. 

By  varying  the  initial  conductivity,  we  found  the  optimal  tracking 
force  for  a  variety  of  beam/channel  parameters.  In  Table  1,  the  conduc¬ 
tivity  which  produced  the  maximum  tracking  force  for  a  conductivity  chan¬ 
nel  is  given.  The  nominal  parameters  are: 

e  =  0.75  cm, 

rc  =  1  cm, 

ab  =  0.5  cm, 

IQ  =  10  kA, 

and , 

1/0  =  15  cm. 

In  Table  1,  the  beam  and  channel  are  always  nominal  unless  specified  other- 
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The  electric  and  magnetic  components  of  the  total 
tracking  force  are  plotted  for  a  10-kA  beam  with 
15-ns  rise  time  in  a  I-cm  density  channel. 


TABLE  1.  OPTIMIZATION  OF  CHANNEL  TRACKING  IN  CONDUCTIVITY  CHANNEL 


Varied  Parameter 

and  Its  Value 

Optimal  on 
(cm'1) 

Maximum  Tracking 
Force  (gauss) 

Nominal 

4.7 

24.8 

ab  =  0.25  cm 

4.4 

46.3 

ab  =  1.0  cm 

5.0 

4.5 

IQ  =  1  kA 

3.0 

3.2 

IQ  =  5  kA 

3.5 

13.2 

r  =  2.0  cm 

1.6 

17.5 

c 

r  =  3.0  cm 

1.0 

11.3 

c 

1/0  =  30  cm 

4.0 

15.3 

1/0  =  60  cm 

3.5 

9.0 

e  =  0.25  cm 

3.6 

00 

• 

'-■J 

e  =  0.50  cm 

4.2 

17.0 

c  =  1.50  cm 

5.0 

36.6 

9  A 


The  values  for  optimal  Og  and  F^  are  surprisingly  high.  For  nominal 
parameters,  a  tracking  force  of  25  gauss  for  4.7-cm"*  conductivity  was  cal¬ 
culated.  As  predicted  in  the  analytic  model  derived  in  the  second  section, 

2 

F^  increases  approximately  as  e Iq e/rc  .  The  optimal  Og  changes  particu¬ 
larly  when  rc  is  varied.  This  phenomenon  is  not  fully  modeled  by  the 
analytic  expression,  since  the  beam  was  restricted  to  be  entirely  within 
the  channel  and  conductivity  generation  was  not  considered.  We  see  the 
optimal  Og  decreases  by  a  factor  of  5  as  rc  increases  from  1  to  3  cm. 
Another  result  not  predicted  by  the  model  is  the  ab  dependence.  The 
tracking  force  increases  an  order  of  magnitude  as  ab  decreases  from  1  to 
1/4  cm,  probably  because  a  greater  portion  of  the  Bennett-profi le  beam  lies 
within  the  channel  for  a  smaller  Bennett  radius.  With  the  maximum  tracking 
force  over  40  gauss,  these  results  are  much  more  optimistic  than  those 
previously  documented. 


In  Table  2,  the  IPROP  results  are  summarized  for  the  density  channel 
case.  Here,  the  optimized  parameter  is  the  displacement.  As  seen  by  Hui 
and  Lampe7,  the  beam  will  not  track  until  the  channel  is  approximately 
1/2  rc  off  axis.  The  optimal  e  for  nominal  parameters  is  1.5  cm.  This 
value  varies  little  unless  rc  or  ab  are  changed.  We  find  the  optimal  dis¬ 
placement  scales  as. 


'opt 


rr  + 

c  b 


(32) 


This  scaling  is  caused  by  the  peak  electric  field  occurring  near  the  edge 

of  the  beam.  Since  the  avalanche-produced  conductivity  grows  as  a  function 
2  2 

of  E  /n  ,  the  tracking  occurs  only  when  the  beam  is  well  off  axis  allowing 
one  beam  edge  to  lie  near  the  channel  center.  There  is  a  similar  scaling 
of  Ft  with  ab ,  IQ  and  rise  time  as  in  the  conductivity  case  with  tracking 
about  half  as  great.  For  optimal  displacement,  IPROP  calculates  a  15  gauss 
tracking  force. 


Finally,  we  examine  the  density/conductivity  channel  in  Table  3.  The 
results  are  similar  to  those  of  the  conductivity  channel  with  one 
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TABLE  2.  OPTIMIZATION  OF  TRACKING  IN  DENSITY  CHANNEL. 


Varied  Parameter 
and  Its  Value 

Optimal  e 
(cm) 

Maximum  Tracking 
Force  (gauss) 

Nominal 

1.5 

14.8 

a.  =  0.25  cm 
b 

1.2 

21.0 

afe  =  1.0  cm 

2.0 

4.3 

IQ  -  1  kA 

1.5 

0.39 

IQ  »  5  kA 

1.5 

CO 

• 

r*— 

rc  =  2  cm 

2.5 

14.8 

rc  =  3  cm 

4.6 

5.9 

1/6  =  30  cm 

1.5 

9.1 

1/6  «  60  cm 

1.5 

3.5 

TABLE  3.  OPTIMIZATION  OF 

Varied  Parameter 

and  Its  Value 

TRACKING  IN 

Optimal  an 
(cm’1) 

DENSITY/CONDUCTIVITY  CHANNEL. 

Maximum  Tracking 
Force  (gauss) 

Nominal 

4.5 

20.6 

ab  =  0.25  cm 

3.2 

33.6 

a.  =  1.0  cm 
b 

5.1 

5.6 

I0  =  1  kA 

2.9 

3.1 

IQ  =  5  kA 

4.3 

13.5 

rc  =  1.5  cm 

1.6 

13.0 

rc  =  2.0  cm 

0.84 

8.2 

1/8  =  30  cm 

4.3 

14.2 

1/8  =  60  cm 

3.4 

8.5 

c  =  1.5  cm 

4.6 

34.4 

e  =  2.0  cm 

4.8 

37.1 

28 


3 


I 


qualification.  At  e  =  0.75  cm,  the  detracking  effects  of  the  density  chan¬ 


nel  near  the  axis  produce  smaller  F^.  We  see  that  when  the  displacement  is 


increased  to  1.5  cm,  the  tracking  force  nearly  catches  up  to  the  conduc¬ 
tivity  case.  Clearly,  at  these  high  conductivities  the  detracking  effects 
of  the  depressed  density  in  the  c r  jination  channel  make  it  less  desirable 
than  the  pure  conductivity  channel. 


B. 


TRUMPET  BEAM 


The  trumpet-shaped  beam  is  a  more  realistic  representation  of  a 
laboratory  beam  and  must  be  given  some  attention  here.  We  will  not  go  into 
as  much  detail  as  in  the  pencil  case,  since  the  results  are  qualitatively 
similar.  As  one  would  expect  from  the  beam  radius  dependence  of  a  pencil 
beam,  the  tracking  for  a  trumpet  beam  is  much  smaller.  The  trumpet  shape 
used  in  IPROP  scales  the  Bennett  radius  to  x  as, 


afa(x)  =  0.5  cm  +  c^  exp(-x/c.,) 


(33) 


where  c^  =  6  cm  and  =  8  cm.  The  current  variation  is  still  described  by 
Eq.  27. 


In  Figure  9,  the  tracking  force  of  a  trumpet  beam  in  a  3-cm  channel  is 


compared  with  that  of  a  pencil  beam  with  a.  =  2  cm  for  a  channel  with  a 
-1  D 

0.2-cm  conductivity.  The  2-cm  radius  is  the  average  trumpet  beam  radius 


during  its  tracking  period.  The  maximum  tracking  for  the  trumpet  beam 
occurs  roughly  at  the  same  position  in  the  beam  as  the  pencil  beam  with 
half  the  force.  The  initial  force  for  the  pencil  is  much  greater,  since 
the  trumpet  is  over  3  times  wider.  The  magnetic  detracking  of  the  trumpet 
beam  is  much  stronger  further  back  in  the  beam,  causing  a  quicker  decline 
than  the  pencil  beam  at  x  =  13  cm.  Despite  these  differences,  the  tracking 
behavior  with  respect  to  channel  offset,  radius,  peak  beam  current  and  rise 
time  is  similar  to  the  wide  pencil  beam. 
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Since  its  effective  radius  is  much  greater,  the  trumpet  beam  needs  a 
correspondingly  greater  channel  radius  to  cause  significant  tracking.  In 

Table  4,  the  optimal  Og  and  tracking  force  are  listed  for  =  1  to  8  cm. 

For  the  same  r^  used  in  the  pencil  case  (1  cm),  the  tracking  force  is 

negligible.  However,  a  peak  force  of  0.78  gauss  is  computed  for  re  =  4  cm. 

As  before,  the  optimal  Og  declines  with  increasing  r^  indicative  of  mag¬ 
netic  detracking  effects. 

The  results  for  trumpet  beams  in  density  and  combination  channels  show 
slightly  lower  tracking  than  the  conductivity  channel  except  at  large  dis¬ 
placements  as  seen  in  the  pencil  beam  case.  Since  the  effective  Bennett 
radius  is  much  greater  in  the  trumpet  case,  in  a  density  channel,  detrack¬ 
ing  effects  are  felt  for  larger  displacements.  The  optimal  displacement  is 
correspondingly  larger  at  5  cm  in  a  3-cm  radius  channel.  For  this  offset, 
the  maximum  tracking  force  is  calculated  at  5.6  gauss  in  the  density  chan¬ 
nel  (see  Figure  10)  and  4.8  gauss  in  the  combination  channel  .  In  the 
density  channel,  the  magnetic  and  electric  components  of  the  tracking  force 
cancel  until  x  =  10  cm.  Beyond  this  point,  the  magnetic  forces  dominate 
until  the  electric  component  finally  becomes  tracking  at  x  =  18  cm. 


TABLE  4. 

TRUMPET  BEAM  TRACKING 

IN  CONDUCTIVITY  CHANNEL 

r 

c 

(cm) 

Optimal  on 
(cm'1) 

Maximum  Tracking 

Force  (gauss) 

1 

8.0 

7  *  10'5 

2 

2.8 

0.052 

3 

1.2 

0.51 

4 

0.78 

0.78 

8 

0.70 

0.70 
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FORCE  (GAUSS) 


As  seen  in  Figure  10,  a  trumpet  beam  experiences  tracking  well  into 
the  beam  for  large  displacements  («  1.5  rc).  For  all  three  types  of  chan¬ 
nels,  no  net  detracking  was  calculated  up  to  t  =  20  cm.  This  behavior  is 
not  nearly  as  pronounced  in  the  pencil  simulations.  Since  the  beam¬ 
generated  conductivity  is  more  diluted  in  the  trumpet-beam  case,  detracking 
axial  return  currents  do  not  grow  as  fast.  This  effect  is  greatest  at 
large  di spl acements  since  the  initial  conductivity  (avalanche-induced  o  in 
the  density  channel  case)  near  the  beam  is  small. 


CHAPTER  4 
DISCUSSION 


F 


Analytic  models  and  IPROP  simulations  have  given  us  insight  into  the 
physical  mechanisms  of  channel  tracking.  Electrostatic  tracking  is  a 
secondary  effect  in  a  beam  with  short  rise  times  (<4  ns)  and  Oq  <  5  cm  ^ . 
This  was  seen  to  be  true  particularly  in  the  case  of  a  density  channel, 
where  electrostatic  forces  are  generally  detracking.  Although  the  magnetic 
tracking  must  overcome  the  electric  detracking,  we  still  calculate  very 
high  forces  on  the  order  of  25  gauss  for  nominal  parameters.  When  the  more 
physically  realistic  trumpet  beams  were  considered,  6-gauss  tracking  forces 
were  still  attainable  at  large  displacements. 

The  magnetic  component  of  the  tracking  force  is  driven  by  J  and  J 

r  0 

channel  currents  and  by  the  dipole  axial  displacement  current.  The  track¬ 
ing  currents  are  produced  by  the  growing  transverse  electric  fields.  These 
currents  in  turn  drive  magnetic  fields  which  attract  the  beam  to  the  chan¬ 
nel.  The  dipole  axial  displacement  current  is  caused  by  the  beam  seeing  an 
increasingly  asymmetric  channel -electron  profile  as  t  increases  for  T  < 

1/a.  After  a  fast  initial  rise  near  the  beam  head,  these  fields  decay 
roughly  as  exp(-ax),  as  predicted  in  the  analytic  model. 

The  main  parameter  which  determines  whether  a  beam  will  track  in 
either  a  conductivity  or  density  channel  is  the  amount  of  volume-integrated 
conductivity  in  the  channel.  This  is  evident  from  the  lower  optimal  aQ 
calculated  for  larger  radii  channels.  This  conductivity  coupled  with  the 
changing  beam  electric  field  produces  a  channel  return  current.  This 
oppositely  flowing  current  will  attempt  to  detrack  an  off-center  beam. 

For  beam-weighted  a  above  1-2  cm'1  in  the  nominal  case,  both  the  analytic 
expression  and  IPROP  predict  magnetic  detracking. 

The  detracking  axial  return  currents  grow  as  a  builds  up  in  the  chan¬ 
nel.  This  current  growth  is  the  result  of  the  initial  conductivity  and 
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i  impact  ionization  of  the  beam.  Close  to  the  beam  head,  E  ^  couples  with 

the  asymmetric  conductivity  that  was  present  initially  or  produced  via 
electron  avalanche  (in  density  channel)  to  produce  dipole  axial  return  cur¬ 
rents.  Also,  the  beam-generated  conductivity  couples  with  E  The 

•  detracking  force  produced  at  a  radius  r  from  the  beam  (r+  is  the  point 

towards  the  channel  axis  and  r"  is  away  from  the  channel)  grows  as, 

Fdt(r)  -  ( a+  -  a’)Ez(0)(r)  +  2a1Ez(1)(r)  (34) 

where  is  the  beam-symmetric  conductivity  produced  in  impact  ionization 
at  r  from  the  beam  axis.  a+  and  a~  are  the  asymmetric  conductivity  at  r+ 
and  r”,  respectively.  For  small  t,  E  ^  scales  with  dlz/dr  and  Ez^ 
scales  with  agd^/d-r.  Since  the  magnitude  of  the  asymmetric  conductivity, 
with  a  correction  for  beam  and  channel  radii,  roughly  determines  (o+  -  a"), 
both  terms  in  Eq.  34  increase  with  increasing  <Jq.  We  therefore  expect 
smaller-og  channels  and  longer  rise  time  beams  (c^  grows  slower)  to  track 
1 onger. 

As  seen  in  Tables  1-3,  there  is  a  strong  increase  in  the  tracking 
force  as  the  beam  radius  declines  for  a  given  beam  current.  This  a^  depen¬ 
dence  is  caused  by  a  greater  fraction  of  the  beam  propagating  in  the  chan¬ 
nel  for  the  smaller  beam  radius.  Initially,  a  thin  beam  experiences 
enhanced  magnetic  tracking  due  to  the  large  localized  electric  fields. 
However,  for  a  given  oq,  this  beam  will  be  repulsed  sooner,  since  will 
grow  quickly  close  to  the  beam  axis.  For  this  reason,  smaller-ag  channels 
are  optimal  for  smaller-a^  beams.  Since  a  greater  percentage  of  the  cur¬ 
rent  lies  in  the  channel  where  it  is  useful  for  tracking,  a  thin  beam  in  a 
channel  with  small  eg  experiences  maximum  channel  tracking. 

The  trumpet-shaped  beams  behave  similarly  to  the  pencil  beam,  if  we 
consider  the  beam  to  have  a  larger  effective  radius.  The  trumpet  beam 
therefore  requires  a  larger  channel  radius  for  good  tracking.  In  the  case 
we  considered  with  an  effective  radius  of  2  cm,  the  ideal  channel  had 
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rc  =  4  cm.  This  is  similar  to  the  scaling  of  the  pencil  beam  which  has  a 

1.2  cm  ideal  channel  with  a,  =  0.5  cm.  This  effect  does  tend  to  lower  the 

b  -2 

maximum  tracking  due  to  the  r scaling  of  the  tracking  force. 
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CHAPTER  5 
CONCLUSIONS 


Maxwell's  complete  equations  have  yielded  very  optimistic  estimates 
to"  channel  tracking.  Both  analytic  and  numerical  results  show  that  the 
large  magnitude  of  the  tracking  is  the  result  of  the  magnetic  forces  pro¬ 
duced  near  the  beam  head.  These  forces  rise  quickly  before  the  electro¬ 
static  forces  have  had  a  chance  to  build  up.  The  main  detracking  force, 
occurring  at  roughly  t  =  5  cm  for  the  nominal  beam/channel  parameters,  is 
also  magnetic,  resulting  from  channel  return  currents.  The  tracking  force 
is  as  high  as  6  gauss  for  a  trumpet  beam  and  up  to  50  gauss  for  a  pencil 
beam. 

The  effect  of  most  beam/channel  parameters,  in  the  IPROP  studies,  is 

similar  to  that  predicted  in  the  analytic  model  derived  in  Section  2.  The 

2 

peak  tracking  force  scales  as  eIQ B/rc  for  ab,  e  <_  in  a  Gaussian  channel 
with  initial  conductivity.  The  tracking  forces  seen  in  IPROP  simulations 
are  severely  degraded  for  a^  _>  rc ,  with  maximum  tracking  occurs  for  e  - 
ab  +  V  In  the  densit*  and  combination  channels,  the  asymmetric  conduc¬ 
tivity  produced  through  electron  avalanche  causes  detracking  for  small  e 
(<  r^/2).  However,  for  larger  e,  the  tracking  force  scales  similarly  to 
the  conductivity  channel. 

Because  of  the  large  transverse  dipole  magnetic  fields,  a  larger 
initial  channel  conductivity  produces  optimal  tracking  than  previously 
thought.  For  a  pencil  beam  in  a  1-cm  radius  channel,  the  maximum  tracking 
occurs  with  <jQ  =  4.7  cm”^,  an  order  of  magnitude  higher  than  the  0.2  value 
suggested  by  Hui  and  Lampe7 .  The  optimal  conductivity  declines  with 
increasing  channel  radius,  since  channel  return  currents  can  occur  in  a 
larger  volume,  causing  detracking.  The  optimal  conductivity,  therefore, 

9 

scales  as  1/r  . 
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In  general,  a  trumpet  beam  behaves  as  a  pencil  beam  with  a  larger 
effective  radius.  The  maximum  tracking  is  achieved  in  a  correspondingly 
larger  channel  giving  roughly  an  order  of  magnitude  less  force  than  a  0.5- 
cm  radius  pencil  beam.  One  significant  deviation  from  pencil  behavior  is 
the  longer  duration  of  the  force  in  a  trumpet  beam.  Because  the  small  cur¬ 
rent  density  near  the  beam  head  produces  a  diluted  beam-generated  conduc¬ 
tivity,  the  tracking  force  persists  farther  back  into  the  trumpet  beam 
particularly  for  e  _>  r  . 

The  results  presented  here  are  preliminary,  since  they  were  calculated 
assuming  a  slug  beam.  In  this  case,  the  tracking  forces  were  easy  to 
distinguish.  In  a  dynamic  beam,  oscillations,  particularly  from  the  hose 
instability,  can  obscure  these  forces.  In  the  near  future,  we  will  con¬ 
sider  the  more  physical  representation  of  an  electron  beam  using  IPROP  with 
dynamic  particles  and  hopefully  better  estimate  the  tracking  force.  The 
work  presented  in  this  paper  does  suggest  we  will  find  significantly  higher 
forces  than  previously  expected. 
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